parpool(6);
pathdir='mesh/';
i=1;j=2;
path1=num2str(i);path2=num2str(j);
tic;
% [vertices1, faces1 ]= obj__read( [pathdir,'mesh',path1,'.obj'] );
% [vertices2, faces2 ]= obj__read( [pathdir,'mesh',path2,'.obj'] );
[vertices1]=readobj_vertex([pathdir,'mesh',path1,'.obj']);
[vertices2]=readobj_vertex([pathdir,'mesh',path2,'.obj']);
toc
P=vertices1';Q=vertices2';
tic;
[d, D] = HausdorffDist(P,Q,0);
toc

tic;
[d, D] = HausdorffDist_par(P,Q,1);
toc

% delete(gcp('nocreate'))